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METHODE POUR DEFORMER GRADUELLEMENT DES SIMULATIONS SEQUENTIELLES D'UN MILIEU 
HETEROGENE TEL QU'UNE ZONE SOUTERRAINE. 

© - Methode pour deformer graduellement les represen- 
tations ou realisations, generees par simulation sequentiel- 
le, d'un modele stochastique non necessairement gaussien 
d'une grandeur physique z dans un milieu heterogene 
maille, afin de les ajuster a un ensemble de donnees relati- 
ves a la structure ou I'etat du milieu qui sont collectees par 
des mesures et observations prealables. 

- Elle comporte essentiellement I'application d'un algo- 
rithme de deformation graduelle d'un modele stochastique 
a un vecteur gaussien a N variables mutuellement indepen- 
dantes, qui est relie a un vecteur uniforme a N variables uni- 
formes mutuellement independantes par la fonction de 
repartition gaussienne de faf;on a definir des realisations du 
vecteur uniforme, et I'utilisation de ces realisations pour ge- 
nerer des representations de cette grandeur physique z, 
que Ton cale par rapport aux donnees. 

- Applications par exemple pour visualiser la configura- 
tion statistique d'une grandeur: la permeabilite d'un gise- 
ment souterrain, la pollution atmospherique, etc. 
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La presence invention a pour objet une methode pour deformer 
graduellement des representations ou realisations, gene'res par simulation 
sequentielle, d'un modele stochastique non necessairement gaussien d'un 
milieu heterogene, en se basant sur un algorithme de deformation graduelle de 
modeles stochastiques gaussiens. 

La methode selon 1* invention trouve des applications dans la 
moderation de zones souterraines ou il s'agit de generer des representations 
montrant comment est distribute une certaine grandeur physique dans une 
zone du sous-sol (la permeabilite z par exemple), qui soient compatibles au 
mieux avec des donnees observes ou mesurees : donnees geologiques, 
enregistrements sismiques, mesures obtenues dans des puits notamment 
mesures des variations au cours du temps de la pression et du debit de fluides 
issus d'un gisement etc. 

Etat de la technique 

Par la demande de brevet FR 98/09.018, on connait une methode pour 
deformer graduellement un modele stochastique (de type gaussien ou 
apparente) d'un milieu heterogene tel qu'une zone souterraine, contraint par 
une ensemble de parametres relatifs a la structure du milieu. Cette methode 
comporte le tirage d'un nombre p (p=2 par exemple) de realisations (ou 
representations) independantes du modele ou d une partie au moins du modele 
choisi du milieu, parmi l ensemble de toutes les realisations possibles et une 
ou plusieurs etapes iteratives de deformation graduelle du modele en 
effecruant une ou plusieurs combinaisons lineaires successives de p 
realisations initiales independantes enue elles puis des realisations composites 
successivement obtenues eventuellement avec de nouveaux tirages, etc.. les 
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coefficients de cette combinaison etant tels que la somme de leurs carres est 
egale a 1 . 

Les modeles gaussiens ou apparentes conviennent pour modeliser des 
champs de grandeurs continues et, de ce fait, sont mal adaptes a la 
5 modelisation de zones traversees par des reseaux de fractures ou des systemes 
de chenaux, par exemple. 

On connait par ailleurs les algoridimes dits de simulation sequentielle, 
les plus largement utilises parmi les algorithmes de simulation geostatistique. 
Bien qu'ils conviennent particulierement pour la simulation des modeles 
to gaussiens, il n'impliquent pas en principe une limitation a ce type de modele. 

On forme une representation geostatistique d'une zone souterraine par 
exemple en la subdivisant par un reseau a N mailles et en determinant un 
vecteur aleatoire a N dimensions Z = (Z,, Z 2 ,. Z N ) correspondant au mieux a 
des mesures ou observations obtenues sur la zone. Comme le montre par 

15 exemple Johnson, M.E. ; in « Multivariate Statistical Simulation » ; Wiley & 
Sons, New York, 1987, cette approche reduit le probleme de la creation d'un 
vecteur a N dimensions a une serie de N problemes mono-dimensionnels. Un 
tel vecteur aleatoire n'est ni necessairement multi-gaussien ni stationnaire. La 
simulation sequentielle de Z implique d'abord la definition d'un ordre selon 

20 lequel les N elements (Zj, Z 2 ..... Z N ) du vecteur Z sont generes Tun apres 
I* autre. Sans rien perdre de generalite, on suppose que les N elements de Z 
sont generes en sequence depuis Zi jusqu'aZ.vj. Pour tirer une valeur de chaque 
element Zi. (i = 1, X), il faut realiser les operations suivantes : 

a) construire la distribution de Z^onditionnee par (Zj. Z:.. Z t .\) 
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F c (z,) = P(Z i <z i /Z l ,Z 1 Z : .,) ;et 

b) tirer une valeur de Zj de la distribution F c (z;). 

Dans la pratique de la geostatistique, la simulation sequentielle est 
utilisee frequemment pour engendrer des vecteurs multi-gaussiens et des 
5 vecteurs d'indicateur non-gaussiens. La tache principale de la simulation 
sequentielle est la determination des distributions conditionnelles F c (z.) (i = I, 
.... N). Des algorithmes et logiciels pour estimer ces distributions sont decrits 
par exemple dans : 

- Deutsch, C.V. et al, « GSLIB (Geostatistical Software Library) and User's 
io Guide » ; Oxford University Press, New York, Oxford 1 992. 

En ce qui concerne le tirage des valeurs a partir de la distribution F c ( Zj ). 
il existe aussi un vaste ensemble d'algorithmes connus. 

On considere la methode de distribution inverse par laquelle on obtient 
une realisation de Z, : = F^(u,), ou Ui est tire d'une distribution uniforme 
15 entre 0 et 1. II s'ensuit qu'une realisation du vecteur Z correspond a une 
realisation du vecteur U dont les elements U„ L%, ....,U N , sont mutuellement 
independants et uniformement distribues entre 0 et 1 

On s'apercoit qu'une simulation sequentielle est une operation S qui 

transforme un vecteur uniforme U= (U„ U : ,U N ) en un vecteur structure Z 

20 =(Z,.Z : ,....Z N ) : 

Z = S(U) (|) 

Le probleme que represente la contrainte d un vecteur Z a differents 
types de donnees peut etre resolu en contraignant les distributions 
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conditionnelles F c (zO (i = 1, N) et/ou le vecteur uniforme U= (U|, U:, 
....,Un). 

Des travaux recents sur Falgorithme sequentiel ont ete concentres sur 
T amelioration de I'estimation des distributions conditionnelles F c (zO par 
integration de donnees geologiques et de donnees sismiques. 

On peut citer par exemple le papier de Zhu, H. et al ; « Formatting and 
Integrating Soft Data : Stochastic Imaging via the Markov-Bayes Algorithm » 
in Soares, A., Ed. Geostatistics Troia 92, vol.1 : Kluwer Acad. Publ., 
Dordrecht, The Netherlands, p. 1-12, 1993. 

Cependant, cette approche ne peut etre etendue a T integration de 
donnees non lineaires comme les pressions issues des tests de puits et les 
historiques de production, sauf h imposer une severe linearisation. En outre, 
comme toute combinaison de vecteurs uniformes U ne donne pas un vecteur 
uniforme, la methode de deformation graduelle d'un modele stochastique 
developpde dans le brevet dejh cite, ne peut etre directement appliquee dans le 
cadre de la technique sequentielle rappelee ci-dessus. 

Aussi, la methode selon 1* invention permet de rendre compatible les 
deux approches i.e. d'etendre le formal isme developpe dans le brevet 
precedemment cite, a la deformation graduelle de realisations, generees par 
simulation sequentielle, d'un modele stochastique non necessairement 
gaussien. 

Definition de la methode 

La methode permet de deformer graduellement une representation ou 
realisation, generee par simulation sequentielle, d'un modele stochastique non 
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necessairement gaussien d'une grandeur physique z dans un milieu heterogene 
tel qu'une zone souterraine, afin de la contraindre a un ensemble de donnees 
collectees dans le milieu par des mesures et observations preambles, relatives a 
son etat ou sa structure. 

Elle est caracterisee en ce qu'elle comporte V application d un 
algorithme de deformation graduelle d'un modele stochastique a un vecteur 
gaussien (Y) ay ant un nombre N de variables mutuellement independantes qui 
est relie a un vecteur uniforme (U) a N variables uniformes mutuellement 
independantes par la fonction de repartition gaussienne (G). de facon a definir 
une chame de realisations u(t) du vecteur (U), et Iutilisation de ces 
realisations u(t) pour generer des realisations z(t) de cette grandeur physique 
que Ton cale par rapport aux donnees (non lineaires). 

Suivant un premier mode de mise en oeuvre, on definit la chafne de 
realisations u(t) du vecteur uniforme (U) a partir d'une combinaison lineaire 
de realisations du vecteur gaussien (Y) dont les coefficients de combinaison 
sont tels que la somme de leurs Carre's est egale a 1'unite. 

Suivant un autre mode de mise en ceuvre. on realise une deformation 
graduelle d'un nombre n de parties du modele represents du milieu 
heterogene tout en preservant la continuite entre ces n parties du modele. par 
une subdivision du vecteur uniforme (U) en un nombre n de sous-vecteurs 
mutuellement independants. 

D'autres caracteristiques et avamages de la methode selon Imvention 
apparaitront a la lecture de la description ci-apres dun exemple non limitatif 
de realisation, en se referant aux dessins annexes ou : 
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la Fig.l montre la couche mediane (Tune realisation d % un modele de facies 
servant de reference, generee par simulation sequentielle dMndicatrices ; 

la Fig. 2 montre la variation en fonction du temps de la pression obtenue 
dans un test de puits du modele de reference ; 

les Fig.3A h 3E montrent respectivement cinq realisations initiales de la 
couche mediane d'une zone reservoir, contraintes seulement par les facias le 
long du puit ; 

les Fig. 4A a 4E montrent respectivement pour ces cinq realisations, les 
courbes des pressions au puits dans le modele de reference comparees a 
celles correspondant aux modeles initiaux ; 

les Fig.5A & 5E montrent respectivement cinq realisations de la couche 
mediane du modele de facies conditionnees aux facies le long du puits et 
calees par rapport a la courbe de pression obtenue par test de puits ; 

les Fig. 6A a 6E montrent respectivement pour les cinq realisations, les 
courbes des pressions au puits dans le modele de reference comparees & 
celles correspondant aux modeles cales ; 

les Fig.7A h 7E montrent respectivement comment les fonctions objectif 
correspondant respectivement a ces cinq exemples, varient avec le nombre 
d" iterations ; 

les Fig. 8A a 8E montrent les deformations graduelles operees par un 
changement de coefficient d'anisotropie, sur un modele a trois facies genere 
par simulation sequentielle d'indicatrices; et 
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- les Fig. 9A a 9E montrent les deformations graduelles locales d'un modele a 
trois fades, generees par simulation sequentielle d'indicatrices. 

Description detaillee de la methode 

On considere une zone d'etude que Ton subdivise par une grille a N 
mailles et Ton cherche a construire des realisations ou representations d'un 
modele stochastique d'une certaine grandeur physique z representant par 
exemple la permeabilite des formations dans la zone. Le modele recherche doit 
s'ajuster a des donnees obtenues par des mesures ou observations en un certain 
nombre de points, et notamment a des donnees non lineaires. 



Le calage d'un modele stochastique a des donnees non lineaires peut 
etre formulee comme un probleme d' optimisation. On designe par f bs = (f, ohs , 
f ^obs ^ fj obs ^obs^ Ie vecteur des jo^^j non ii n ^ a i res observees ou mesurees 

dans le milieu etudie (la zone resen'oir), et par f = (f, , f : , f 3 , f p ), j e 

vecteur correspondant des reponses du modele stochastique de la permeabilite 
id Z = (Zi, Z;,..., Z N ). Le probleme de contraindre le modele stochastique Z par 
les observations consiste a generer une realisation z de Z qui reduit a une 
valeur assez faible une fonction objectif que Ton definit comme la somme des 
erreurs quadratiques ponderees des reponses du modele par rapport aux 
observations ou mesures dans la zone-reserv oir, soit : 



20 o=^5>(/,-rv 



oil co, represente le poids attribue a la reponse f,. Les fonctions f, (i=l. 2, 
p) ainsi que la fonction objectif O sont des fonctions du vecteur Z. On se 
trouve devant un probleme d'optimisation de dimension N. 
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Pour etendre le formalisme developpe' dans le brevet deja cite, a la 
deformation graduelle de realisations generees par simulation sequentielle non 
necessairement gaussienne, on va partir d'un vecteur gaussien a N variables Y, 
avec i = 1, 2, N, mutuellement independantes, de moyenne nulle et de 
variance egale a 1' unite, et on va definir N variables uniformes mutuellement 
independantes Ui, U ; , U 3 ,...Un par : 

Uj = G(Yj) V i = 1,2, .... N 

ou G designe la fonction de repartition gaussienne standardisee. 

Ceci pose, on va appliquer Talgorithme de deformation graduelle 

developpe dans un cadre gaussien, au vecteur gaussien Y = (Y,, Y 2 , , Y N ). 

pour construire tine chaine continue de realisations du vecteur uniforme U = 

(Ui, U2, Un). Etant donne deux realisations independantes y a et y b de Y, on 

definit la chaine de realisations u(t) du vecteur U que I on obtient selon la 
relation : 

m(0 = G(y a cos/ + y b sin/) (2) 

Pour chaque t, u(t) est une realisation du vecteur U. Ensuite. par 

echantillonnage de la disu-ibution conditionnelle F.(z,)(i=l, 2, N) en 

utilisant les elements du vecteur u(t). on obtient un vecteur z(t) qui est, pour 
chaque t, une realisation du vecteur aleatoire Z. En consequence, on peut 
ajuster le parametre t comrae dans le cas gaussien de facon a caler z(t) a des 
donnees non lineaires. Cette procedure est iteree jusqu'a atteindre un calage 
satisfaisant. 
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Calage d'un modele de fades a des donnees de pression obtenues par 
des tests de puits 

Pour illustrer I' application de la methode d' optimisation stochastique 
ci-dessus definie, on va chercher a caler un modele stochastique de reservoir a 
des donnees de pressions obtenues par des tests de puits. La construction du 
modele de reservoir est derivee d'une formation petroliere reelle qui comporte 
trois types de fades : deux fades reservoir de bonne qualite (fades 1 et 2) et 
un type de fades reservoir de tres mauvaise qualite (fades 3). La table I 
definit les proprietes petrophysiques des trois facies : 
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Pour representer la distribution specifique de facies de la formation 
petrolifere, on ge'nere d'abord une realisation binaire pour representer le facies 
3 et son complement. Ensuite, dans la partie complementaire du facies 3, on 
genere une autre realisation binaire independante de la premiere pour 
representer les facies 1 et 2. La formation est discretisee par un maillage 
regulier de 60x59x15 blocs de taille !5mxl5mx!.5m. Un modele de 
variogramme exponentiel est utilise pour estimer les distributions 
conditionnelles. La direction d' anisotropic principale est diagonale par rapport 
au maillage. Les portees du variogramme du facies 3 selon les trois directions 
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d' anisotropic sont respectivement de 300m, 80m et 3m. Les portees du 
variogramme des fades I et 2 selon les trois directions d'anisotropie sont 
respectivement de 150m, 40m et 1.5m. Les proportions des facies [, 2, 3 sont 
respectivement de 6%, 16% et 78%. 

On a realise un test de puits au moyen d'un simulateur de test de puits a 
differences finies, de la mani£re decrite par : 

Blanc, G. et al : « Building Geostatistical Models Constrained by Dynamic 
Data - A Posteriori Constraints » in SPE 35478, Proc. NPF/SPE European 3D 
Reservoir Modelling Conference, Stavenger, Norway, 1996. 

On voit sur la Fig. 1 la couche mediane d'une realisation qui sert de 
modele de reference pour cette validation. La section du puits qui a ete 
perforee traverse horizontalement la couche mediane du modele de reservoir le 
long de Taxe x. Le diametre du puits est de 7.85cm, la capacite du puits est 
nulle et les facteurs parietaux (skin) des facies 1, 2 et 3 sont respectivement 0, 
3 et 50. Le test de puits synthetique dure 240 jours avec un debit constant de 5 
m3/jour de fagon a investiguer presque la totalite du champ petrolifere. La 
Fig. 2 montre la variation de pression en fonction du temps. 

Le but etait de construire des realisations du modele de facies contraint 
par les facias rencontre le long du puits et par la courbe de pression obtenue 
lors du test de puits. La fonction objectif est definie comme la somme des 
differences quadratiques entre les reponses en pression du modele de reference 
et les reponses en pression de la realisation. Du fait que le comportement 
dynamique du modele de reservoir est principalement controle par le contraste 
entre les facies reservoir de bonne et de mauvaise qualites, la realisation 
binaire utilisee pour gdn<frer les facies I et 2 a ete fixee d'abord et seule la 
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realisation binaire utilisee pour ge'nerer le fades 3 a ete deformee pour le 
calage des donne'es de pression. 

Les reponses en pression resultant des tests de puits pour les cinq 
realisations des Fig.3A a 3E, sont differentes de celle du modele de reference, 
conune le montrent les Fig. 4A 5 4E. Partant respectivement de ces 5 
realisations independantes, en utilisant la methode de calage iteratif ci-dessus, 
on obtient apres quelques iterations, cinq realisations calees (Fig. 5A a 5E) 
pour lesquelles les courbes de pression correspondantes sont tout a fait en 
accord avec celles du modele de reference, conune le montrent les Fig.6A a 
6E. 

Deformation graduelle par rapport aux parametres structuraux 

Dans de nombreux cas, on ne dispose pas de donnees suffisantes pour 
en deduire les parametres structuraux du modele stochastique : moyenne, 
variance, fonction de covariance, etc. Ces parametres structuraux sont souvent 
donnes en termes d' interfiles ou de distributions a priori. Si leurs valcurs sont 
fausses, il est vain de rechercher une realisation calee sur des donnees non 
lineaires. Aussi, il est essentiel pour les applications de pouvoir proceder a une 
deformation graduelle d'une realisation avec modification simultanee de 
nombres aleatoires et des parametres structuraux. L'algorithme de simulation 
sequentielle definie par I'equation (1) fournit cette possibilite de changer 
simultanement ou separement Toperateur structurel S et le vecteur uniforme U. 
Les Fig.8A-8E montrent les deformations graduelles que Ton obtient pour une 
realisation fixe du vecteur uniforme V quand on change le coefficient 
d' anisotropic 
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Deformation graduelle locale ou regionalisee 

Quand les observations sont reparties sur des zones differentes d'une 
formation etudiee, un calage qui utiliserait une deformation globale risquerait 
d'etre inefficace car r amelioration de 1' accord que Ton obtiendrait sur une 
zone pourrait le deteriorer sur une autre zone. II est done bien preferable 
d'appliquer des deformations graduelles zone par zone. On considere une 

subdivision du vecteur U en un certain nombre n de sous-vecteurs U 1 . U" 

U n mutuellement independants ce qui rend possible de realiser 
individuellement leur deformation graduelle. En appliquant separement a 

chacun des sous-vecteurs U 1 , U 2 , , U n , ralgorithme de deformation 

graduelle, on obtient une fonction de dimension n, du vecteur uniforme U : 

G^cost, ♦ ^sinf,)" 
G(Y; cos t 2 + Y;s\nt 2 ) 

G(Y;cost n + Y»iint n \ 

ou KJ et Yl pour tout i =1, 2, n, sont des sous-vecteurs gaussiens 
independants. Pour un ensemble de realisations de K,' et Y*. on resout un 
probleme d* optimisation de n parametres tj, t : , t n pour obtenir une 
realisation qui ameliore, sinon maintient, le calage des donnees. Cette 
procedure peut etre iteree jusqu'i atteindre un calage satisfaisant. 

En realisant des deformations locales graduelles, on ameliore de faqon 
significative la vitesse de calage dans tous les cas ou les mesures ou 
observations sont reparties dans differentes zones du milieu. 



U z (t z ) 
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On voit bien Teffet de cette deformation locale graduelle sur le modele 
a trois fades, des Fig. 9A i 9E ou seule la partie inferieure gauche delimitee 
est affectee. 

La methode selon 1' invention peut etre facilement generalisee h la 
deformation graduelle d.'une representation ou realisation d un modele 
stochastique quelconque, car generer une realisation d'un tel modele 
stochastique se ramdne toujours a la generation des nombres uni formes. 
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REVINDICATIONS 



\ 

1) Methode pour deformer graduellement une representation ou 
5 realisation, gener^e par simulation sequentielle/d'un modele stochastique non 

necessairement gaussien d'une grandeur physique z dans un milieu heterogene 
tel qu'une zone souterraine, afin de la contraindre a un ensemble de donnees 
collectees dans le milieu par des mesures et observations prealables, relatives i 
son etat ou sa structure, caracterisee en ce qu'elle comporte F application d'un 

io algorithme de deformation graduelle d'un modele stochastique a un vecteur 
gaussien (Y) a N variables mutuellement independantes qui est relie a un 
vecteur uniforme U a N variables uniformes mutuellement independantes par 
la fonction de repartition gaussienne (G) ? de fagon i construire une chame de 
realisations u(t) du vecteur U, et T utilisation de ces realisations u(t) pour 

15 generer des realisations z(t) de cette grandeur physique, que Ton cale par 
rapport aux donnees. 

2) Methode selon la revindication 1, caracterisee en ce que Ton definit 
une chaine de realisations u(t) du vecteur (U) a partir d'une combinaison 
lineaire de realisations du vecteur gaussien (Y) dont les coefficients de 

20 combinaison sont tels que la somme de leurs carres est egale a Tunite. 

3) VIethode selon la revendication 1 ou 2, caracterisee en ce qu'elle 
comporte une deformation graduelle du modele representatif du milieu 
heterogene, simultanement par rapport aux parametres structuraux et aux 
nombres aleatoires. 
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4) Methode selon la revendication 1 ou 2, caracterisee en ce qu'elle 
comporte une deformation graduelle separee d'un nombre n de parlies du 
modele representatif du milieu heterogene tout en preservant la continuite 
entre ces n parties du modele, par une subdivision du vecteur uniforme en n 
sous-vecteurs mutuellement independants. 
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Fig. 1 Fig.2 
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Fig.6A Fig.6B Fig.6C Fig.6D 
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